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We have evaluated the Drell-Yan cross section in nuclei paying spe- 
cial attention to the meson cloud contribution from pion and p-meson, 
for which an accurate calculation using the meson nuclear spectral func- 
tions is used. Similarly, the nucleonic contribution is evaluated in terms 
^ ' of a relativistic nucleon spectral function. Fair agreement with exper- 

Ch ' iment is found for different nuclei and the results show a sizeable con- 

tribution from the renormalized meson cloud. In order to reproduce 
i^ ' the experiment a novel element is introduced, consisting of a gradual 

'j_j ■ energy loss of the incoming proton in its pass through the nucleus which 

5t ! produces a strong A dependence at xi large. 



PACS: 24.85. +p, 13.85.Qk, 14.40.Ag, 25.40.Ep 



The Drell-Yan process (DY) in nuclei in which a quark (antiquark) from 
a proton beam and an antiquark (quark) from the constituents of the nuclear 
target fuse to produce a /i+^u" pair, has been advocated as a complement of 
deep inelastic scattering (DIS) in nuclei in order to investigate the reasons why 
the quarks in a nuclear medium (even if they are integrated into other effective 
degrees of freedom, like nucleons and mesons) behave differently than when 
they are in isolated nucleons |l|, ^, |^, ||, |^, ^ |^. Comprehensive reviews of this 
process and its relation to the EMC effect can be seen in |^, ^, |TU|, |TT|. 



Since in the DY process one needs an antiquark for the reaction, there was 
the hope that, compared to the DIS, the nuclear DY would be more sensitive to 
modifications of the pion cloud in the nucleus, one of the mechanisms originally 
suggested to explain the EMC effect ||, |, |10|, |ll|. 



The calculations in ref. show that the region of large xi and small X2 
in the Drell-Yan process is enhanced due to medium modifications of the pion 
cloud, while the region of small xi would simply reproduce the same results 
as the EMC effect. In ref. [^ the pionic excess leads to a sharp increase of the 
DY ratio for large Xi and X2 < 0.5. Similarly an increase of that ratio with 
increasing values of Xi ■ X2 is predicted in due to mesonic effects. 

The experimental data, however, does not show the expected enhancement 



12| and recent papers try to justify the small role of the pion cloud in this 



process and consequently in the EMC effect [0. On the other hand, the 
mesonic contribution in the EMC effect was revived recently in [Q, where a 
new evaluation was carried out in terms of the pion spectral function, avoiding 
the static approximation (use of pion excess number) and demanding that the 
pion propagator satisfies rigorous analytical properties. 

The mesonic effects of ref. |T^ showed up as an enhancement of the ratio 



F2A/F2D in the region 0.1 < a; < 0.5, in qualitative agreement with previous 
findings. However, the contribution from the renormalized pion cloud was 
found smaller than claimed in the past [^ |T0[. At the same, time the renor- 



malized p-meson cloud in the medium was shown to be as relevant as the one 
of the pion, and even more important in some kinematical regions. Together 
with the use of an accurate nucleon spectral function, including relativistic 
effects, the analysis of |T^ showed a good agreement with experiment in the 



EMC region for a wide range of nuclei. 

In the present work we would like to extend the ideas of ref. [n| to the 



Drell-Yan process and establish the role of the meson cloud in that reaction. 
The results show an important role of the meson cloud in the DY process but, 
unlike in earlier works based upon the meson cloud, we can now find agreement 
with experiment. In ref. [0 the effects of the pion cloud in the DY process 



are also reanalyzed demanding similar requirements as those of ref. \Ti\. The 
effects of the p-meson cloud are not considered there. 

There is another important novelty in our work with respect to former 
works in the DY process. Since the protons of the beam are strongly interacting 
particles, they will have strong collisions with the nucleons much more often 
than they will eventually produce the ^~^fi~ pair. We shall take this into 



consideration here. We shall work outside the shadowing region for the variable 
X2 in order to avoid the typical coherent phenomena of this region. Outside 
this shadowing region the different parts of the nuclear volume will contribute 
incoherently to the cross section. The incoming proton collides strongly with 
the nucleons, but the proton keeps travelling and there is no loss of flux. 
However, in any of these collisions the proton of the beam will lose a certain 
amount of energy. This will change the variable Xi and consequently the 
contribution to the cross section. Let us see how this is implemented: 
The cross section for pN -^ /i"'"/i~X is given by 

(fa{pN -^ /iV~^) = n 2 I]ga[ga(3:i)ga(a;2) 

QQ a 

+ qa{Xi)qa{x2)] dXidX2 (1) 

where the sum is over the flavours of the quarks and Xi,X2 refer to the beam 
and target nucleons respectively. The variables Xi,X2 indicate the fraction 
of the momentum carried by the quark (antiquark) of the beam and target 
nucleons which fuse to create the virtual photon that leads to the fi'^fx^ pair 
of momentum q. In an invariant form we have 

Xi = 2 ; X2 = 2 ; sxiX2 = q (2) 

s s 

where pi,P2 are the fourmomenta of the beam and target nucleons and s = 

{pi +^2)^- In the frame where the target nucleon is at rest we can write: 

where M is the nucleon mass and Ei the energy of the incoming nucleon. The 
axis 3 is chosen along the direction of the beam. 

On the other hand we shall have a distribution of j92;P2 i^ the nucleus given 
by the nuclear spectral function. This will give rise to a variable X2n which in 
the Bjorken limit can be written in terms of the static X2 variable as 

X2N = 2 = X2-0 o (4) 

sn m - P2 



which is the same relationship of xn to a; in D IS p4 |. 

The amount of energy that the nucleon of the beam loses in one collision 
with the nucleons in the nucleus is difficult to quantize. On the one hand some 
energy is transferred to the target and many new particles can be created 
at the high energies which we are discussing, Ei ~ 800 GeV. On the other 
hand some times the nucleon can get excited to some resonant state in the 
strong collision and the resonance may behave like a nucleon with respect 
to the electromagnetic /i^/i~ process. Furthermore, given the scale of time 
in which the process occurs, the /i"*"//" production could take place before the 
asymptotic final state from the strong collision materializes. All this tells about 



the difficulties in determining the equivalent energy loss for one collision. For 
this reason we do not attempt to evaluate this magnitude, but use the same 
DY experiment to fix it. We assume that in each collision, occurring with a 
probability aNNpdl, we lose a fraction /3 of the energy. This fraction does not 
need to be the same for all energies, but we will assume it to be constant in the 
energy range where we move. We take the same function for computations in 
different nuclei. The value of (3 is taken such as to reproduce one experimental 
point in one nucleus, which will be described later on. Thus we have: 



^ = -ctnnpPE, 

Ei{r) = Eiin exp[-/5(TAriv/_!^p(6,z')c?^'] 



(5) 



xi{r) = xi exp[PaNN Hoc Pib, z')dz'] 

where a^^ is the A^A^ total cross section, a^N = 40 mb, p{r) the nuclear 
density and b the impact parameter. We have assumed that the nucleon loses 
energy but keeps moving in the forward direction, in the spirit of the eikonal 
approximation. In proton elastic collisions with the nucleon {aei — (Ttoi/^) the 
amount of energy lost is negligible because the cross section is very forward 
peaked [l^. So, strictly speaking we should use (Tj„ instead of aNN, but since 
what matters is the product /3 ■ a and (3 is fitted to the data, we can keep using 
the formalism of eq. (5). 

With this prescription and the formalism of ref. [l^ we obtain 



rfV(^) {j)A -^ fi+fi-X) A-Ko^K 



X 



dxidx2 
d^p M 



A d\Y.^l[qa{xi{r)) 



{2TifE{p)J- 
d^p M r^^ 



9g2 



dp^Sh{p^,p)qa{x2N) + qa{xi{r)) 



dp''SK{p\p)qa{x2N)]e{x2N)e{l~X2N)9{l-Xi{r)) (6) 



(27r)3^(p) 
where Sh{jf,p) is the hole spectral function for the nucleon in the nucleus. 



For the pion cloud contribution we have, following ref. |T4 

rfV^'^) {pA -^ fi+fi-X) Ana'^K 



dxidx2 



9g2 



-6) f d^r 



IIe^[g„(xi(f)) 

a 

+qaixi{r)) 



d'^p 



9{pyimD^{q)2Mqa{x2^] 



d^p 



(2vr)4- 

e{p°)SImD^{q)2Mq,{x2.)] 



(2vr)^ 

9{x2n - X2) 9(1 - X2n) 9(1 - Xi(f)) 

The p contribution is obtained from eq. (7) changing D-^{q) 
X2p and the factor (— 6) by (— 12). 



(7) 



-Dp(g), X2„ 



The magnitude SD{p) is given, following ref. [14|, by 



SD{p) = D{p) - Doip) 



dD 

dp 



■ P 



p=0 



with D{p) standing for the it ot p propagator in the medium and -Do(p) the 
corresponding free one. Furthermore 



M 



X2tt — X2 



P2-P2 



(9) 



P2 referring now to the pion momenta. The expression of X2n compared to 
the one of X2n in eq. (4) has an apparent minus sign which comes because 
the variable p for pions has opposite direction as for nucleons for convenience. 
The variable X2p is given by eq. (9) but now p2 stands for the p momentum. 
The extra factor 2 appearing for the p meson contribution in the cross sec- 
tion, with respect to the expression for pions, accounts for the two transverse 
polarizations of the p which couple to the nucleons in our approach |T^ . 

In eq. (6) the q,q with arguments Xi{f) refer to the proton beam and 
<ia{,X2N) refer to the nucleon target (average over proton and neutron). 

In eq. (7), instead, the g, g with argument X2tt{x2p) refer to the quark 
distributions of the pion (p meson, this latter one is taken the same as for 



pions ||I6|, |171). The quark distributions of nucleons and pions are taken from 
refs. 



1^ 



and [0 respectively. For the meson propagator and nucleon spectral 
function we use the same input as in the study of the EMC effect in ref. ||14|| . 

A small technical change is made here. The form factor accompanying 
the coupling of the mesons to the ph, Ah components is taken here of the 
monopole type as in [|1^], with the same values of A, but we take it static 
(dependent on three momentum only) to avoid unrealistic contributions from 
regions in the integration where one is close to the poles of the form factor. 
We should note that these form factors are constructed in the study of the 
NN interaction where g^ < |20]. While this does not practically modify the 
results in DIS, it produces some reduction of the results here. Particularly, 
the meson contribution to the DY in one nucleon, which can be obtained by 
substituting 5D{p) by the last term in eq. (8), is more strongly reduced and 
reaches reasonable numbers. 

Next we show some results. We have evaluated the results for R = 2daA/AdaD 
corresponding to the experiment of [O . We have evaluated the results in two 
different ways. First we take a value of X2 and evaluate a weighed cross sec- 
tion multiplying the cross sections at different values of Xi by the experimental 
acceptances ||21|. As indicated in |]12|, we include in the average only the re- 
gions 4 GeV < \^CP < 9 GeV and \/Q^ > 11 GeV, in order to avoid the 
regions of quarkonium resonances. Another calculation is done by assuming 
xi— X2 = 0.26 (the average value in [|l^) in order to approximately extrapolate 
the results to the unmeasured region. We do not evaluate da for the deuteron 
in our local density approach and hence we divide by the average of the cross 



section on the proton and the neutron. For values of X2 < 0.6 this has errors 
of less than 2 % as found in the study of DIS p2|, p3 . 



We fix the fraction of energy loss in order to obtain one point of the spec- 
trum of ^^Fe, at X2 = 0.15. The resulting fraction is (3 = 0.035, or 3.5% energy 
loss. This is a small quantity but it plays a role as we shall see. 

In figs. 1,2,3 we show the results for ^^C, '^''Ca and ^^Fe. The experimental 
results are roughly reproduced in the three nuclei. We notice the relative 
importance of using the weighed cross section instead of the one using the 
average value of xi — X2, although the differences are smaller than 10%. The 
latter cross sections show a small A dependence which is different than the 
one of the weighed cross sections. The ratios obtained from the weighed cross 
sections are of the order of unity, as approximately shown by the experiment, 
although the ^^C data seems to be larger than unity around X2 — 0.25. The 
results using the xi — X2 average also show a ratio bigger than unity for ^^C 
in this region. 

It is interesting to show the effects of the meson cloud. In figs. 1,2,3 we 
have separated the nucleonic contribution from the one of the pion and p- 
meson clouds. As we can see, the contribution of the pion cloud is a bit larger 
than the p one. The role of the energy loss can be seen in fig. 1. At small 
values of X2, taking xi — X2 = 0.26, it reduces the cross section in about 10 %. 
At larger values of X2 the reduction is much bigger. This stronger reduction is 
easily understood since, due to the energy loss, Xi increases from its original 
value X2 + 0.26 and eventually becomes bigger than unity where there is no 
strength for the cross section. 

It is interesting to see, however what are our predictions in two limiting 
cases, Xi small and xi large. In fig. 4 we show the results for the three nuclei 
for Xi = 0.01 and Xi = 0.7. For xi = 0.01 we obtain values of R very similar to 
those obtained for the EMC effect, as already noticed in 0. The A dependence 
in this case is very weak. However, for xi = 0.7 the dependence of R on X2 
is quite different for all three nuclei. It is easy to understand this different 
behaviour: for xi very small, the energy degradation increases xi by a certain 
fraction but it is still very small and this energy loss induces small changes 
(smaller than 10 %). On the contrary, when Xi is close to unity, a fractional 
increase of Xi brings it closer to one where there is no strength. This induces 
large reductions in the cross sections, which are larger in heavier nuclei, hence 
the A dependence seen in the cross sections. 

Although our pionic effects are qualitatively similar to those in earlier 
works, the energy loss changes appreciably the spectra with respect to for- 
mer predictions [@, ^. 

The issue of the energy loss in hadronic processes is attracting some atten- 
tion. Early estimates of quark energy losses based on the uncertainty principle 
pointed towards very small energy losses of quarks propagating through the 
nucleus ||2^. More accurate estimates, yet not free of uncertainties as claimed 



by the authors, are done in ES 



These estimates are lower than 3.5% energy loss per collision. More recent 

6 



evaluations [^ distinguish between quarks created in the nuclear medium 
and incoming quarks. For the latter case an energy loss proportional to the 
energy is obtained which can be translated in our language as about 2.5% 
energy loss per collision. On the other hand, the energy loss per collision 
counted asymptotically, is certainly larger than 3.5% in order to reproduce 
the broad energy distributions of the experimental cross sections [^ (see also 



comments in ref. |^ and particularly fig. 8). These experimental facts can be 
reconciled recalling our arguments that in some cases hadronic resonances are 
excited, carrying the energy of the incoming proton, which propagate through 
the nucleus and can produce /i"*"//" similarly as nucleons. However, in pure 
hadronic reactions these resonances would decay outside the nucleus into a 
nucleon and mostly pions and the nucleon would have less energy than the 
incoming one. 

It is also interesting to quote that in a recent paper p9| looking at the 
propagation of J/ip in nuclei a solution was favoured implying both absorp- 
tion of J /if) through the nucleus and energy loss of the beam. This loss was 
equivalent to 3.6% per collision, although some small trade-offs could be made 
between the absorption and the energy loss. 

Given the relevance of the mesonic components and the energy loss in the 
DY process, we should worry about questions of selfconsistency since the pion 
structure functions are determined from analysis of DY processes in which no 
energy loss is assumed. A look at fig. 4 would tell us that for values of Xi 
small, where the energy loss played a small role, one has not much to worry, 
but for large values of Xi such things could be more relevant. This would add 
certain uncertainties to our predictions for xi = 0.7 in fig. 4. 

In any case the study done here clearly shows that the issue of the energy 
loss is an important one and experimental efforts should be devoted to clarify 
it. Fig. 4, even with accepted uncertainties, also shows that the DY process 
at large values of Xi is the relevant place to look at. 

In summary, our results show a sensitivity of R to the meson cloud renor- 
malization in nuclei, more important than in the EMC effect. The values of 
R around unity in our interpretation are not a signal of the lack of mesonic 
effect, but the simultaneous effect of the mesonic cloud and the progressive 
energy degradation of the nucleon beam through the nucleus. The latter has 
the effect of reducing the cross section since xi increases and q{xi) {q{xi)) 
decreases. The present interpretation of the DY nuclear effect has as a con- 
sequence a stronger nuclear dependence than in DIS, where the ratio is prac- 
tically constant for A > 7. Here we have found that for xi close to unity 
the ratio R is rather dependent on A. It would be very interesting to have 
other DY experiments done with a larger range of values of a;i,a;2 and better 
precision in order to be able to test the novel consequences that the present 
interpretation of the process provides. 
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Figure Captions 

• Figure 1: The ratio 2daA/ ^dao for tlie Drell-Yan process on ^^C Ex- 
perimental points from ref. [0. Long-dashed line: nucleonic contribution. 
Dashed-dotted hne: plus pionic contribution. Solid line: plus p-meson contri- 
bution (full calculation). The former curves use the experimental acceptances 
and cuts. Short dashed line: full calculation using xi — X2 = 0.26. Dotted line: 
full calculation with xi — X2 = 0.26 omitting the energy loss. 

• Figure 2: Same as fig. 1 for ^°Ca. 

• Figure 3: Same as fig. 1 for ^^Fe. 

• Figure 4: Three upper curves around small X2'- values of R for fixed value 
Xi = 0.01 for ^^C, ^''Ca, ^^Fe (from up down), as a function of X2- The lower 
curves around small X2'- values of R for fixed value Xi = 0.7 for ^^C, ^°Ca, ^^Fe 
(from up down), as a function of X2. 
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